PROCEEDINGS 

OF SCIENCE 



Lattice simulation of 

ultracold atomic Bose-Fermi mixtures 



Arata Yamamoto 

Quantum Hadron Physics Laboratory, RIKEN Nishina Center, Saitama 351-0198, Japan 



E-mail: arayamamoto@riken.jp 



Bose-Fermi mixtures have been recently realized and invesitigated in ultracold atomic experi- 
ments. We formulate quantum Monte Carlo simulation of Bose-Fermi mixtures on the (3+1)- 
dimensional lattice. As its first application, we analyze the boson-fermion pair correlation and the 
phase diagram of the Bose-Einstein condensation. 



The XXX International Symposium on Lattice Field Theory 
June 24-29, 2012 

Cairns Convention Centre, Cairns, Australia 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Lattice simulation of ultracold atomic Bose-Fermi mixtures 



1. Introduction 

One frontier of ultracold atomic physics is multi-component quantum system, such as Bose- 
Fermi mixtures. The most famous example is a 4 He- 3 He mixture. Recently, many kinds of Bose- 
Fermi mixtures have been investigated in ultracold atomic experiments. In particular, Bose-Fermi 
mixtures can be trapped on a three-dimensional optical lattice. The mixtures of 87 Rb- 40 K [jlj ^[ ^J, 
170 Yb- 173 Yb and 174 Yb- 173 Yb [Q] were examined on the optical lattice. Since the optical lattice is 
a real physical lattice, we can exactly study lattice physics in laboratory. 

On the theoretical side, we can exactly reproduce the lattice physics in quantum Monte Carlo 
simulations. Bose-Fermi mixtures have been studied in quantum Monte Carlo simulations in 1+1 
dimensions [Q, |6[ [7|, ||, ^, [K], 11]. All the previous simulations were done in the world-line formal- 



ism [12] or its extensions [13, 14]. The world-line formalism is an exact scheme in 1+1 dimensions. 
In higher dimensions, however, there is the sign problem which originates from antisymmetric 
property of fermions. For this reason, there was no quantum Monte Carlo study of Bose-Fermi 
mixtures in 3+1 dimensions. 

In this study, we perform quantum Monte Carlo simulation of Bose-Fermi mixtures on the 
(3+l)-dimensional lattice. We adopt the same framework as the lattice QCD simulation. In this 
framework, there is no sign problem on fermions. We use the lattice unit and drop the lattice 
constant throughout this paper. 

2. Formalism 

We consider one-component boson field <J>(x, t) and two-component fermion field ^(x, x) 
and \Pj_(jc, f). The formalism is based on the path integral in terms of the Euclidean action. The 
generating functional is 

Z = J D®*D$D x ¥\D x ¥ t D x ¥lD x ¥± e~ s 

= J D®*D<& det/Y t det7Y ; ^ SB , (2.1) 

This system is similar to the two-flavor QCD, in which the boson field is gluon and the fermion 
fields are u-quark and d-quark. We adopted the hybrid Monte Carlo algorithm, which is frequently 
used in the lattice QCD simulation. 

For the lattice action, we consider the Bose-Fermi Hubbard model. The naive form of the 
action is 

S = S B + S F + S BF (2.2) 
Sb = £ <I>*(x,t){<I>(x,t)-<I>(x,t-1)}-^«b(x,t) 

,V,T 

3 

-t B £ {<&*(*, x)&@+ej, t) + <&*(*, r)<S>(x-ej, t)} + U B n B (x, z){n B (x, z) - 1} (2.3) 

7=1 
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3 

-t Pa T)^ ff (x + ? 7 -,T) T)V a (jt-2j, T)} (2.4) 

Sbf = £ U B pn B (x,r)n Fa (x,r) (2.5) 

X.T,(7 

with n fl (jf, t) = <J>*(x, t)<&(x, t) and n/?cr(jc, t) = *F£(*, t^^x, t). 
We list several remarks on this action: 

• This action has the sign problem on the imaginary-time derivative term of the boson field. 
This is understood in the Fourier transformation as 

t){0(x,t)-0(*,t-1)}=W^^ . (2.6) 

The boson Matsubara frequency is (0 B k = 2knT. To avoid this sign problem, we adopted the 
zero-frequency approximation. In this approximation, nonzero-frequency modes <&(x,k ^ 0) 
are set to zero, and the boson field is independent of imaginary time as <&(x, t) = &(x,k = 
0) = <J>(x). This approximation is a priori justified in high-temperature limit or near the 
critical temperature. In general case, we can check the validity of the approximation by 
adding a few lowest frequency modes, e.g., in the reweighting method. This sign problem is 
characteristic in the non-relativistic theory. There is no sign problem in the relativistic scalar 
theory because the imaginary-time derivative term is bilinear. 

• We do not consider the fermion self-interaction term for the simplicity. It is straightforward 
to treat the fermion self-interaction term by the Hubbard-Stratonovich transformation and 
the auxiliary field. 

• We take the same hopping parameter and the same chemical potential for the two fermions, 
tp = tp^ = tp± and }Xp = jUf-j- = Hfi, because the hybrid Monte Carlo algorithm becomes 
simple. The two fermions are degenerated. 



The imaginary-time derivative is discretized to the backward difference. The chemical poten- 



tial and the interaction are multiplied to the backward hopping term, according to Refs. [ 15 



16]. For example, 



z)[W(x, t) - Vfc z-1)] + {-Hf + U BF n B (x)}y*(x, z)¥{x, t) 
^V*(x,z)[W(x,z) + {-e^+U BF n B (x)}W{x,z-l)} . (2.7) 



• The most time-consuming part of the simulation is the inversion of the fermion matrix. We 
used the BiCGstab solver for the matrix inversion. We performed the Fourier transforma- 
tion to the frequency space for preconditioning, i.e., for accelerating the convergence of the 
solver. After the Fourier transformation, the fermion matrix becomes diagonal in the fre- 
quency space. Because the temporal lattice size N x is much larger than the spatial lattice size 
N s , this preconditioning greatly improves the convergence of the solver. 
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From the above consideration, the final form of the lattice action is 



Sb = N t *£ 



[l-e^)n B (x 



-tB^{^*{x)^{x + ej) + ^*(x)^(x-ej)} + UBn B (x){n B (x) - 1} 
7=1 



(2.8) 



Sf = N T £ [(l-e^e- i<0n )n Fa (x,k) 

x,k,<y 
3 



-tw £ {% % k)Va (x + e h k)+% (x, k)<¥ a (x- ej , k) } 
5bf = A^t £ UBFe~ ia)Fk n B {x)n Fa {x,k) . 

x,k,a 



(2.9) 
(2.10) 



The fermion Matsubara frequency is a>fk = (2k — l)nT. 

We can prove the positivity and the reality of the fermion determinant. In the fermion matrix, 
the complex factor is only e~ l(0fk and other parts gives real eigenvalues. Therefore, the fermion 
determinant is 



detK a = Y\(A + e- i( ° n B) = ]\ [(A + B cos co Fk ) 2 +B 2 sin 2 co Fk ] >0 

k sin(B fl( .>0 



(2.11) 



where A and B are real numbers. 

There are many simulation parameters in this system. We fixed the hopping parameters 
t B =t F = 0.01, the chemical potentials }X B = }X F = 0, the boson self-interaction U B = 0.1, and 
the spatial lattice volume Nf = 10 3 . We varied the boson-fermion interaction U BF = —0.1 to 0.1 
and temperature T = l/vV T = 0.01 to 0.05. 

3. Result 

Because the simulation is in the grand-canonical formalism, number densities are not fixed. 
The number densities are functions of chemical potentials and other parameters. In Fig. p], we 
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Figure 1: Left: the boson number density Right: the fermion number density (n F ). 
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Figure 2: The boson-fermion correlation function (nB(x)np(y)) as a function of the distance R = \x—y\. 
Left: the attractive boson-fermion interaction Ubf = —0.5. Right: the repulsive boson-fermion interaction 
Ubf = 0.5. 



show the boson number density (ng) and the fermion number density (rip). The number densities 
depend on the boson-fermion interaction Ubf ■ The number densities are decreasing functions of 
Ubf because the chemical potentials are replaced as /ig — > /ig — Ubf^f and pLp —> — Ubf^b at 
tree level. Since the free fermion has the particle-hole symmetry, the fermion is half-filling, i.e., 
(rip) = 0.5 at Ubf = 0. On the other hand, the boson does not have the particle-hole symmetry at 
finite Ub- Thus, the regions of Ubf > and Ubf < are not symmetric. 

In this setup, we measured the boson-fermion pair correlation (nB(x)np(y)). At x = y, this 
quantity is called the pair occupancy. The physical interpretation of the pair occupancy is a proba- 
bility to find boson-fermion pairs in a single lattice site. In Fig. ^ we plot the boson-fermion pair 
correlation as a function of the distance R = \x—y\. The temperature is T = 0.05. In the attractive 
case Ubf < 0, the pair correlation is enhanced at R = 0, and thus the formation of boson-fermion 
pairs is favored. In the repulsive case Ubf > 0, the pair correlation is reduced at R = 0, and thus 
the bosons and the fermions tend to separate. 

Next we analyze the Bose-Einstein condensation. The Bose-Einstein condensation is identified 
from the long-range behavior of the boson propagator (x)<&(y)), which is so-called the off- 
diagonal long-range order. As shown in the left panel of Fig. ||, the boson propagator drops to zero 
at high temperature. At low temperature, zero-momentum mode of the boson field appears and the 
boson propagator has a nonzero expectation value in R — > 00. This expectation value corresponds 
to the Bose-Einstein condensation density. In the right panel of Fig. ||, we draw the condensate 
fraction 

(n B ) (&(x)<t>(y)) R=0 ' 1 " ' 

Despite the fixed boson self-interaction Ub, the phase transition temperature is changed by the 
boson-fermion interaction Ubf ■ We see that the fermion-induced interaction affects the Bose- 
Einstein condensation. This is similar to the confinement-deconfinement phase transition in QCD. 
Color confinement is a gluon phenomenon, but its phase transition temperature is changed by the 
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Figure 3: Left: the boson propagator (<t>' (x)<t>(y)} as a function of the distance R = \x — y\ at Ubf = 0. 
Right: the phase diagram of the Bose-Einstein condensation. 

dynamical quark effect. The current simulation was performed in a fixed spatial volume. For 
precisely determine the phase boundary, we need to perform finite size scaling. 

4. Summary 

We have formulated the quantum Monte Carlo simulation of Bose-Fermi mixtures on the 
(3+l)-dimensional lattice. We have calculated the boson-fermion pair correlation and the phase 
diagram of the Bose-Einstein condensation. Physics in 3+1 dimensions differs from physics in 1+1 
dimensions. For example, the Bose-Einstein condensation cannot be observed in 1+1 dimensions. 
We should note that these are experimental observables in ultracold atomic physics. The pair oc- 
cupancy was experimentally measured on an optical lattice The fermion-induced effect on the 
Bose-Einstein condensation was observed in the interference pattern on an optical lattice [Jxj] - 
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